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It has recently been shown that second-order corrections to the background distance-redshift 
relation can build up significantly at large redshifts, due to an aggregation of gravitational lensing 
events. This shifts the expectation value of the distance to the CMB by 1%. In this paper we show 
that this shift is already properly accounted for in standard CMB analyses. We clarify the role that 
the area distance to the CMB plays in the presence of second-order lensing corrections. 


I. INTRODUCTION 


Second-order perturbations in the standard model are important for precision cosmology. They carry with them 
a difficulty in interpretation, because they do not average to zero, and consequently give a (small) change to the 
background. The most important are potential changes to observables: are these ever significant, and can they lead 
to biases in cosmological parameter estimation? This question is the essence of the problem of backreaction. Analyses 
in recent years have shown that ‘backreaction’ probably cannot be the origin of dark energy, but it has to be taken 
into account for precision cosmology as it could bias cosmological parameters at the level of up to several percent, see 
e.g., 0" 


Clearly, like all other cosmological quantities the area distance is fluctuating and the power spectrum of its first 
order perturbations has been determined (Toi. [Til]. Recently it was shown that second-order lensing shifts the area 
distance (also called angular diameter distance) to the CMB. The expectation value (d*) of the area distance to the 
CMB is larger than the background value do* = X */(1 + z *)> [l2j: 


(d*) — do* 


! + §<*?> 


( 1 ) 


where K\ is the first-order magnification at redshift z*. At large distances the shift in the mean area distance builds 
up proportional to the comoving volume, approximately given by: 
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( 2 ) 


where x* denotes the comoving distance to z*. This is a 1% increase for canonical cosmological parameters. (Here 
and throughout we only consider the dominant contributions - there are smaller effects such as time delay and 
ISW contributions in every quantity we calculate, see discussion in Section IIII Bl Within this approximation the 
relative change in the area distance and in the luminosity distance is the same, therefore we shall not systematically 
distinguish between them but just use the term ’distance’ which can always be replaced by ’area distance’.) Naively, 
we can associate a perturbative change in the distance (d) = do(l + A) with a change in the spectrum to 

C e =(^) C^ e = (l-2A)C {1 _ Ay , (3) 

where Cf, is the spectrum in the background [13]. For A > 0 as in eq. m this gives a shift in the peaks to larger 
i. Assuming that the shifted distance is measured by the CMB, this increases the estimation of Hg by 5%, while 
decreasing fl m by 10% fl2j . 
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However, CMB observations are not directly sensitive to the expectation value of the distance but rather to its 
average over directions. This is closer in essence to how we extract a model from observations - by averaging over the 
sky. The angular power spectrum of the CMB, C(_, can be approximated by averaging the temperature anisotropy over 
many patches of size 7r/£ in different directions of the sky. In a companion paper [141 ] . we show that at second order 
in perturbation theory averaging over directions followed by an expectation value is not equivalent to the expectation 
value (followed by an average over directions). An alternative, but more relevant, mean distance is therefore 


^ d* ^ ^ ^ ^ dfl no d* (n D ) ^ — do* 


1 - 2 ^> 


( 4 ) 


Lensing therefore also generates a decrease in the observed angular mean distance to the last scattering surface. It is 
remarkable that although the expectation value of the distance along a single line of sight is increased, the total effect 
on the observers sky serves to cancel this effect, and bring the mean distance last scattering surface closer. Interpreted 
as in eq. ©. this shifts the peaks to lower i. 

In this paper, we show that eq. © is too simplistic to describe the change in the spectrum when the magnification 
matrix is affected by a combination of shear and convergence. We derive an approximate expression for the lensed 
angular power spectrum as a function of the convergence and the shear, including second-order contributions. We 
show that the convergence at second order does not contribute to the averaged Ci since it can be written as a total 
divergence. Hence lensing changes the averaged Ce via the square of first-order terms only, and produces a shift to 
lower £. 

We compare our derivation with standard CMB derivations in terms of the bending angle a. Standard expressions 
take into account a variety of lensing effects from terms of the form af ~ O(nf), where oq is the first-order bending 
angle. We extend this calculation to the second-order bending angle ol 2 and show that this term does not contribute 
to the averaged C?. We conclude therefore that no second-order terms have been neglected in standard CMB analyses 
and that lensing cannot be responsible for the tension between local and CMB measurements of Ho as suggested could 
be the case in [12 ] . 

The remainder of the paper is organised as follow: in section [© we present the standard derivation of the lensed Ce 
in terms of the deflection angle a, extending it up to second order. In section HTT1 we derive an approximate expression 
for the lensed angular power spectrum as a function of the convergence and the shear at second order. We conclude 
in section Hvl 


II. THE CONTRIBUTION OF OL TO THE LENSED POWER SPECTRUM 

Lensing shifts the position of the points on the last-scattering surface. The lensed temperature observed in direction 
n Q is therefore given by the unlensed temperature in direction n 


T(n 0 ) = T(n) = T(n 0 + a) ~ T(n 0 ) + a • VT(n Q ) + a J d t djT(n a ) H- (5) 

where a is the bending angle. In the current calculations of the lensed power spectra, only the first-order contribution 
to the deflection angle a is taken into account, both in the second and third term on the right-hand side of eq. © 
0. Here we include the effect of a at second order a = aq + 0:2/2. 

The second-order 02 generates a new contribution to the lensed power spectrum given by the correlation of 


6 2 T = a 2 -VT(n 0 ), 


( 6 ) 


with the unlensed temperature. The deflection angle 02 is derived in Appendix [Cj It reads 

a%(x*,n 0 )=8 [ dx [ dx'— —— —d b d a ^>{x, n 0 )d b ^(x', n 0 ), (7) 

Jo Jo X 

where a = 6 : <p denotes the two components of a transverse to the photon direction n Q . Here d a = e*9,;, with e a 
orthogonal to n D (see appendix C) and 'L denotes the Weyl potential defined in appendix A. 

In Appendix E we show explicitly that the contributions of 8 2 T to the two-point correlation 

(82 atm at m' ) ( 8 ) 

exactly vanishes if we consider ol 2 and T to be uncorrelated. This is a good approximation, since most of the deflection 
angle is generated at z < 100, whereas the CMB anisotropies stem from the last scattering surface (apart from the 
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ISW term which is relevant only on large scales, dominated by cosmic variance). We derive this result both in the 
full-sky and in the flat-sky approximation. 

Alternatively, the result can be understood by noting that 

(<5 2 T(n 0 )T(n')) = <a 2 (n 0 ))<VT(n 0 )T(n' 0 )). (9) 

Since both a 2 and VT(n 0 )T(n' 0 ) are vectors on the sphere, their expectation values have to vanish in order not to 
break statistical isotropy. 

This shows that if we are interested in the effect of lensing on the CMB up to second order in perturbation theory, 
it is sufficient to include the first-order bending angle squared. In section Hill we show how this result translates in 
terms of the convergence and the shear. 

Going to higher order we see that (( 5 iT(n 0 )( 52 T(ng)) contains an odd number of correlators and therefore it vanishes 
for Gaussian initial perturbations. The first non-vanishing higher order terms are therefore (S 2 T(n 0 )S 2 T(n' 0 )) and 
(5±T(n 0 )53T(n' 0 )), which have been calculated in [l 8 |. As is to be expected, the effect on the temperature anisotropies 
of these terms is very small. Nevertheless, for percent accuracy in the polarization, they have to be taken into account; 
in particular, their effect on converting E- into -B-polarization seems to be considerable. However, these results seem 
to contradict a lensing calculation of post-Born terms which obtains much smaller corrections [19j |: therefore, this 
clearly requires more work. 


III. THE CONTRIBUTIONS OF K AND 7 TO THE LENSED POWER SPECTRUM 


We derive now an approximate expression for the lensed angular power spectrum in terms of the shear and the 
convergence. We follow the approach of mm, keeping all relevant terms up to second order in perturbation theory. 
In this approach, the convergence k and the shear 7 = 7 W + are assumed to be constant over the patch of the 
sky in which we measure the temperature fluctuations. This approximation is well motivated to compute the effect 
of a large-scale lensing mode on a small-scale temperature fluctuation. It clearly breaks down when one considers 
fluctuations of the lensing potential at the same scales as the temperature fluctuations. 

Since we are mainly interested in small angular scales, we work in the flat-sky approximation throughout. More 
details on the flat-sky approximation are found in Appendix ID 21 We measure fluctuations in the temperature at two 
observed positions x„ and y D in the sky. These positions are deflected by lensing so that their true positions are x 
and y. If the separation between the points is small we have 


: - y = A(x 0 - y D ), 


( 10 ) 


where A is the magnification matrix 

A = 


1 — k — — q( 2 ) 

—'y ( 2 ) 1 — K + 


= (1 - k )1 + 7 . 


Perturbations in the distance are given by the determinant of A via 

j- = v / det A — 1 k i| 7| 2 , 
do ^ 

where 

| 7| 2 = | dety| = (7 (1) ) 2 + (7 (2) ) 2 ■ 
The lensed CMB power spectrum, C(£), is then given by 


*/f'\\ — f f fyo-iteo+u'yo 


(T(£)T*(£')} = 


2n 


2n 


(T(x 0 )T( y 0 )> = 


f gPxo [ d 2 y p _ lgXn+ityo 


( 11 ) 


( 12 ) 


(13) 


27r 


27r 


(T(x)T(y)} 


f rf2x ° f ^ .-Ox„+il'y, 


27r 
d 2 £ 1 
27r 


27r 


I ^L e «i(x-y ) C (h) = J ^ j flA e -^o+r^yo I ^i e iUA(x 0 -y o)(7(£i) 


C{h)d{£ - A£i)S(£ - £') = - £'), 


(14) 
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where in the second equality we have used that the lensed temperature at position x„ is given by the unlensed 
temperature at position x. For a change in the distance that arises from both convergence and shear then, the correct 
shift in the spectrum is given by 


C(£)=(Jfj C(\A~H\) 


(15) 


and not simply by eq. m, which only applies for a diagonal magnification matrix. 
Keeping terms up to second order in the convergence and the shear we have 


\A-H\ = 


1 + k — i ^ f £ + k 2 + ^|7| 2 — 2 k £ 71 — 71) 2 


t = $t. 


(16) 


The lensed angular power spectrum depends therefore not only on t but also on the direction of £. This breaking of 

isotropy is generated by the terms proportional to £ ')£. 

Defining D{£) = £ 2 C(£), eq. (fl5l) becomes 


D{£) = 


1 


;D(p£) 


detA f3 2 

1 + 2fji- 2 | 7 | 2 + A(£ 7 if + 2k ivy £\ D{£) + 


-£ T i£ + J + \\i \ 2 £D\£) 


+ 


1 9 1 /oT 2 ~T k 

-k 2 + -(£j£) -k£^£ 


£ z D"{£) . 


We are interested in the isotropic part of D(£), averaged over all directions of £. Using that 


If -T ~ 

— dd e £ 7-6 = 0 


and — 
27T 


7 ^) 2 = ^|7| 2 , 


where 9i denotes the direction of £, we find 1 


m= d{£) 


1 , 


K + K + —17| 


£D\t) + 


;l7p 


£ 2 D"(£). 


(17) 


(18) 


(19) 


So far the calculation has been done for a small patch of the sky, where A was assumed to be constant. In practice 
the observed D(£) are effectively averaged over the sky. The expectation value of the averaged spectrum is therefore 
given by 


D(£) = D{t) + 


( k ) + ( K 2 ) + i ( |7| 2 


ID\£) + i 


2 ) + ^ ( l7| 2 


£ 2 D"{£) , 


( 20 ) 


where a bar denotes an average over directions n 0 , as defined in eq. ©• The first order convergence vanishes on 
average (Tcj") = 0. The second-order convergence K 2 can potentially contribute to the average. In appendix E we 
calculate explicitly the convergence up to second order in perturbation theory. We show that K 2 can be written as a 
total derivative 


*= f\fV x ’- xKx - x ' ) 

Jo Jo 


xVx* 


-2 


n o,X)^Kx') + $ ^( n O: x) ^( n O! x' 
n 0 , x) #^(n 0 , x') + ^ 2 ^(n 0 , x) ^(n D , X 


( 21 ) 


where = (4> + ijj)/2 and the transverse operators $ and $ are defined in appendix E see also J25|. As a consequence 
the average of vanishes 2 : (7^2 ) = 0. The only contribution to the lensed spectrum comes therefore from the square 


1 Note that the average over the direction of £. can only be performed after D(f3Z) has been expanded. Averaging directly the argument /3 
in eq. m is not consistent. 

2 Note that for second-order terms, the difference between the average over direction and the ensemble average is higher order so that we 
have (7^2) = (^ 2 )- Recall also that we are only considering terms with 4 transverse derivatives. 
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of the first-order convergence and shear. In appendix [A] we show also that the combination n\ — |q 2 | can be written 
as a total divergence, so that on average 


<«i> = (l7i I) • 


( 22 ) 


With this eq. m becomes 


( D{t) ) = D(£) + \{n\)ZD\l) + 1(4)1* D"(t) . (23) 

The square of the first-order convergence affects therefore the lensed power spectrum. The first term in eq. Cl shifts 
the position of the peaks, whereas the second one both smoothes the peaks and shifts them. These two types of 
corrections are however already consistently included in standard CMB analyses, which include terms up to O (ct 2 ), 
where a is the deflection angle (see section ED- The only contribution not included in previous calculations of the 
lensed CMB spectrum is the second-order convergence Ki. However, as we show here its effect exactly vanishes on 
average, so it does not introduce any additional change in the spectra. 

Finally let us note that if we replace eq. Cl by expression ©, assuming that the changes in the spectrum are due 
only to the mean distance ( d) = do(l - (4 }/ 2 ), instead of eq. (E3l) we would obtain 

(W)) = D{i) + \(4)ZD'(t). (24) 

The shift in smaller than the one given by the first term in eq. (121 and the smoothing term is not present. 


A. Shift in the peaks 

We can now calculate the shift in the position of the peaks induced by the square of the convergence. Denoting the 
(observed) peak of the lensed D(£) by i Q and the peak of the unlensed D(£) by £*, we have 


11 


D'(lo) =0=1 + -<«?> D\£ 0 ) + ~(4)£ 0 D"(t 0 ) + -(4)ilD"\l 0 ). 


Expanding the unlensed spectrum around using £ a = £* + S£ we find 


se 

£* 


(4) 


11 + 3£* 


D"'(l*) 


D"(Q 


(25) 


(26) 


The shift can be calculated by approximating the unlensed spectrum by D(£) oc cos 2 (■k£/£*). Inserting this in 
eq. ([21 gives simply 3 


S£ H('n) 
T* ~ 4~ 


(27) 


Lensing therefore shifts the position of the peaks to smaller multipoles, which is consistent with a decrease in the 
observed distance, as seen in eq. ©. The shift is proportional to (4) which reaches a percent at the last-scattering 
surface. However, the present calculation of the lensed spectrum should be taken with precaution since it is valid 
only for large-scale lensing modes (which can be approximated as constant for a fixed £). The small-scale lensing 
modes contribute however significantly to (4)- For those modes the calculation above cannot be trusted and it is 
necessary to account for convolutions between the temperature and the lensing deflection. The shift calculated in 
eq. d 2 ZD therefore significantly over-estimates the true shift induced by first-order terms squared, but the qualitative 
behaviour is sound. 

Finally as mentioned above, in addition to a shift in the position of the peaks, lensing also induces a smoothing of 
the peaks, due to the second term in eq. d23D- This smoothing actually dominates over the displacement term and it 
constitutes the main impact of lensing on the extraction of cosmological parameters. 


3 Note that since D'''(£*) = 0 , the next term in the expansion around /:* becomes relevant and a more accurate expression for Si in this 


i H( K i) 

case is given by -p ~- . 1 


i _ 3 +i) p2 

1 4 ~ D 7 


+.) 


WT 


_hM) 

4 


[l + 3(K^>7r 2 tt]. 
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B. Why transverse derivatives dominate 


In our derivation of the convergence and the deviation vector at second order, we have neglected all contributions 
with less than four (respectively three) transverse derivatives of the gravitational potential (see appendices lAl and ICl). 
This is justified since the gravitational potential remains small on all scales, while its second spatial derivatives can 
become large, since density fluctuations are large on small scales 

H~ 2 A<j) ~ — > 0( 1) on small scales. (28) 

P 

Moreover, time derivatives of the potential can be neglected with respect to spatial derivatives, since cosmological 
perturbations vary very slowly with time. Finally, radial derivatives on the light-cone are also smaller than transverse 
derivatives. We can indeed rewrite radial derivatives along the null geodesic as 

[ ° dXd r (/)= [ ° d\ n o V0 = [ ° dA(n - a)V0 = 0(A O ) - 0(A„) - [ ° d\ [d t (t> + aV0], (29) 

Jx, Jx, J X, Jx, 

where A* is the value of the affine parameter at the source. As argued before, </> and dt<f>, are much smaller than spatial 
derivatives. In addition, a'Vcj) is a second-order perturbation (since a is itself a perturbation) and it is consequently 
smaller than transverse derivatives of <f>. 


IV. CONCLUSIONS 

In this paper we have derived an approximate expression for the lensed angular power spectrum of the CMB in terms 
of the shear and the convergence. We have consistently included all (dominant) lensing terms up to second order in 
perturbation theory. We have shown that the pure second-order contribution proportional to K 2 where k = Ki + K 2/2 
does not affect the lensed power spectrum, since K 2 can be expressed as a total divergence. Also the second-order 
contribution to the shear only appears at third order. Squares of first-order terms, on the other hand induce a small 
shift of the CMB peaks to lower multipoles. This corresponds to a decrease in the observed mean distance to the 
last scattering surface. We argue that this shift is properly accounted for in standard CMB analyses and that it can 
therefore not be responsible for the tension between local and CMB measurements of Hq as suggested could be the 
case in [ 12 ] . 

The present derivation differs in two points with respect to the analysis in fl2l ] . First, flfjf was using eq. ©, which 
is actually modified when lensing shear is present - we have given the correct version in eq. ED. which is valid for 
a non-diagonal ma gnif ication matrix. Second, jl2| computed the expectation value of the distance, whereas as we 
discuss in detail in [l4| , the quantity which is more relevant to the CMB power spectrum is the expectation value of 
the angular average of the distance. This difference is important, since lensing decreases the angular average of the 
distance, whereas it increases the expectation value of the distance. 

In our derivation we have included only the dominant contributions to the convergence and the deflection angle, 
i.e. those with the maximum number of transverse derivatives. However, since these contributions vanish on average, 
one can wonder if the sub-dominant ones will give an additional shift to the distance. Since the convergence is a 
scalar field, it can only have an even number of transverse derivatives. The terms with no transverse derivatives will 
certainly contribute to the mean distance (since they cannot be written with a total divergence). These terms are 
due for example to the Integrated Sachs-Wolfe effect or the Shapiro time-delay. They affect photon propagation and 
change the physical length of the geodesics between us and the last-scattering surface. Therefore it is natural that 
these terms change the mean distance to the CMB. However, their amplitude is about 1CU 5 at first order and we 
expect the second-order terms to be the square of this. 

The terms with two transverse derivatives may or may not vanish on average. Only a full calculation of these terms 
can determine if they can be written with a total divergence or not. Since these terms describe a coupling between 
longitudinal and transverse deflections, it is probable that they will shift the distance. However, also the impact of 
these terms will be very small. Eq. m shows that first-order terms square (with four transverse derivatives) are of 
the percent level, and therefore terms with only two transverse derivatives will be well below the percent, at most of 
the order of 1CU 6 . (The terms with 2 transverse derivatives can schematically be written as a coupling between the 
first-order convergence and longitudinal perturbations, proportional to 'F, 9 t’F or integrated along the photon 
trajectory. The first-order convergence reaches at most a few percent at the last-scattering surface, whereas the 
longitudinal perturbations are of the order of 10~ 5 .) These subdominant contributions to the mean distance to the 
CMB will therefore probably remain undetected for a long time. 

Note added: Some similar issues are discussed in | 2 Cl | which appeared after this work was completed. 
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Appendix A: The lens map to second order 

We derive the lens map to second order, using the geodesic deviation equation. We follow the formalism and notation 
of [ 24 J], where the shear has been computed to second order. Here we also compute the convergence, however, we only 
consider the perturbations with the maximal number of transverse derivatives. For redshifts z > 0.5 these dominate 
the result. We consider scalar perturbations in longitudinal gauge, 

ds 2 = a 2 (rj) [—(l + 2(jJjdr] 2 + (l — 2?/>)dx 2 ] . (Al) 

Photon propagation is conformally invariant, hence we can ignore the scale factor a{rf). We also introduce the Weyl 
potential *F = (<f> + ip)/2. 

The geodesic deviation equation can be rewritten as an evolution equation for the 2x2 magnification matrix T> ab 

d 2 

~rfj2^ >ab = 'R-ac'Dcb , (A2) 

where lZ ac = Ri_ l v pa k' J k p n£n% and A is the affine parameter along the photon geodesic. Here n p , with a = 1,2, are two 
unit vectors orthogonal to k p and to the observer velocity u p . Eq. (ISZD can be rewritten in terms of the conformal 
distance x = V ~ Vo- At first order, after integration by part, the solution becomes 

'D a b(x*)=[ dx(2 - k°)S ab + [ dx(x* ~ X)xKab, (A3) 

Jo Jo 

where x* is the conformal distance to the last scattering surface. This first-order solution can be used to calculate 
the solution at second order. Formally, at second order we have 

Vab{.X*) = X*fiab + [ dx— - -Sab, (A4) 

Jo X 

where the source term S a b is defined as 

«Mx) = -^y n ^ V ^b - |j^r( 2 - k °) 6 ^ - X^T dx'x'Kab ■ (A5) 

At second order we have two types of terms: the second-order source terms integrated along the background trajectory, 
and the first-order source terms integrated along the perturbed trajectory. Expanding these first-order terms around 
the background trajectory 

Si ab{x pert) = «5i a b(x l ) + 6x[djSi a b, (A6) 

and combining them with the second-order terms, we obtain 

v v dk° d.k° f x {dk°\ 

S ab(x) = Jj^yKacVcb - ~ k °) S ab - X~j^ J dxlxj^-ab + 5x 3 djllab - x&x 3 dj J S ab , (A7) 

where now all the integrals are along the background geodesic relating the image to the observer. The last two terms 
are the corrections to the so-called Born approximation 4 . 


4 Note that the last term is missing in eq. (54) of ref. [23 . Since it only contributes to the trace of X it does not affect the calculation 
of the shear presented there. 
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In the following, we will concentrate on the dominant contributions, i.e. those with four transverse derivatives. The 
terms in lZ ab have at most two transverse derivatives, those in Sx l have at most one transverse derivative, and in k° 
there are no transverse derivatives. From this, we see that the only terms with four transverse derivatives in eq. m 
are the first term, when both !Z ac and V cb are taken at first order; and the fourth term, when both 6x° and TZ ab are 


taken at first order. At first order, the dominant contribution to TZ a b reads 

H lab = -2 eieldidj*, (A8) 

where the vectors e a , with a = 9,ip, are orthogonal to the observed direction n D : 

eg = (cos8cosip,cos9simp, —sin0), (A9) 

e v = (— sin 95 , cosv?, 0 ), (A10) 

n c = (sin 6 cos ip, sin 9 sin ip, cos 9). (All) 

Inserting eq. (1A8I) into eq. (1A8I) . we obtain for the first-order dominant contribution to T> ab 

V lab ( x ) =-2e l a e J b [ d\'(x - x')x'didj®(x') ■ ( Al2 ) 

Jo 

Here and in the following we drop the argument n Q in the gravitational potential when there is no ambiguity. With 
this, the first term in eq. (IA71) becomes 

Siabix) = 4 xe l a e J c e^ei did 3 ^(x) [ dx’(x ~ x!)x’9kde^(x') ■ (A13) 

JO 


To calculate the contribution from the fourth term in eq. (IA71) . we need the deviation vector Sx z at first order. The 
calculation of the deviation vector is presented in appendix [Cj Combining eq. (IC'271) with eq. (1A8[) . we find for the 
fourth term in eq. (E3 


S 2ab (x) = 4x 2 e z c di (eieldjde^ix)) ■ / d\'(x ~ x')e k c d k ^(x') ■ 


(A14) 


The dominant contributions to T> ab at second order then read 


£> 2 ab(x*)= 4 / dx I dx'(x* -X)(X~X') x!e l a e 3 c d l d j ^(x)e k c e f b d k d f .^(x') 


(A15) 


+ Xe l a e J b e k c d l d j d k 'S>(x)e%'S’(x') 


where we have neglected derivatives of the vector e z a , which lead to contributions with less transverse derivatives. 

The transverse derivatives eJ a di can be rewritten in terms of the derivatives on the sphere $ and tf). These operators 
depend explicitly on the spin of the field S X to which they are applied (for more details see e.g. appendix B of 0). 
In terms of the angles 9 and ip, they read 


S X = — sin s 9(dg + i esc 05 ¥> )(sin s 8) S X 


S X = — sin s 9(dg — i esc 09 v )(sin s 9) S X . 


(A16) 


Defining the vectors e± = ± ie v , we have the following relations for a scalar field X = X(x) 


_8 t X = - (jX , 

X 


ei diX = -- $X. 
X 


(A17) 


Then, using that x e ±di e ± = cot 9e± we find 


e\e\did^X = 


l X . 


X 


eLeLdidjX = $X . 

X 


(A18) 


and, analogously with = — cotde± — 2n J 0 , we have 


eleldidjX = e^eldidnX = 


1 w 


- —d r X. 
X 


(A19) 
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Since we are only interested in the contributions with the most transverse derivatives, the last term in eq. (IA19I) can 
be neglected. 

We then decompose the magnification matrix in terms of the shear, the convergence and the rotation 


Dab = A 


* 


q( 2 ) + uj 


where each component has a first and second-order part n = k± 
eq. (1A12I) and using that = (e l + e J _ + e l _e J + )/2 , we obtain 


«i(x*) 

7i(x*) = 7i 1} + wP° 



_^,( 2 ) — 0J \ 

1 — K + 7*- 1 - 1 J ’ 

(A20) 

+ K 2/2 and similarly for 7 and to. 

At hrst order, from 

dx' ' 

(A21) 

XX* 


dx ' 

(A22) 


The rotation u> exactly vanishes at first order. 
At second order, using eq. (I A 1511 , we obtain 


«<*.)= r < ix / v *--* 1 '*-* 0 

Jo J o 


x 2 x'x* 


$ 2 T(x)$ T(x')+ $ ^ (x) $ 2 ’k (x') + 2 W (x) $(xO 


+2 $ + 2 $ 2 ^(x) ^(x') 


(A23) 


These terms can be combined and rewritten with a total transverse derivative, as shown in eq. ED- 

Moreover the combination k 2 — |y | 2 can be computed from eqs. (1A21I) and (IA22I) . It can also be written in terms 
of total transverse derivative 


(«i - |7 iD(x*) = o 


[% /V (x, ~ x) ! x r x,) 

/ 0 j o XXX* 


0(0*(x) WxO- *(x')) 


+ 


')- ^(x) ^(x' 


(A24) 


Appendix B: Explicit calculation of (k\) 


The average of the first order convergence squared can be written as 

(k 2 ) = f dx ( dx'— --— (^$(n 0 ,x)W$(n 0 ,x')) ■ (Bl) 

Jo Jo X*X X x ' 

We can expand the gravitational potential in spherical harmonics 

T(n 0 ,x) = AE*' / d 3 kT*(k,x)aVMkx)YUno)Y; m (k ), (B2) 

where C(k) is the primordial curvature perturbation and T^(fc, x) is the transfer function for the gravitational potential. 
Using that 

A / h\ ns ~ 1 

«(k)C*(k')} = (2n) 3 — J 5(k- k'), (B3) 

where A is the primordial amplitude, n s the spectra index and k\ the pivot scale, and applying the derivative operators 
on the spherical harmonics 

n D ) = -t{l + l)Wm(n 0 ), (B4) 


we obtain 


<« 5 ) 


£ 


-1)£ 2 (£ + l ) 2 


' dk 

~k 


' k_ s 


-1 r 


dx 


X* - X 
XX* 


T*(k,x)je( k X) 


(B5) 
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Appendix C: Deviation vector up to second order 


The deviation vector Sx a is by construction the difference between the position on the perturbed photon geodesic 
and the position on the background geodesic, relating the image to the observer 


6x a = x a — Xq = + ^ x 2 ■ 

The position x a can be calculated from the null geodesic equation 

d 2 X a 

dX' 2 + r ^ k ' J ' k =0 - 

Rewriting A in terms of the conformal distance X = do ~’d and using that 

— . 

dX d\ 


we obtain 


d 2 x a _ dk° k a T “ v ki*k v 


dx 2 


d X (fc 0 ) 2 (fc 0 ) 2 


(Cl) 


(C2) 


(C3) 


(C4) 


At zeroth order, T“ = 0, we have 


1. Order zero 


dfcp 

~dX 


dkg 

dx 


= 0 => kn = cst 


We choose k® = 1 and kg = n l 0 . The position on the background geodesic is then 

x o(X*) = X* ' n o > x o(x*)=X*, 


(C5) 


(C6) 


where we have placed the origin at the observer x“(0) = Xq (0) = 0. 


2. First order 


At first order, (IC4l) reads 


Integrating we obtain 


dk a _ r°j-rk v 
dx k° 


r Q 

L l^v^O^O ■ 


k\x*) 

k°(x*) 


<(l + 2V»(X*)) - 2 [ dxd l ^{x )> 

Jo 


l-2<Kx*) + 2 [ dxd x ^(x), 
Jo 


(C7) 


(C8) 

(C9) 


where we have neglected the perturbations at the observer. Here d x denotes a partial derivative with respect to %. 
We then integrate one more time to find the position x a (x*) : 


gives 


dx a k a 
dx k° 


x l (x*) = x o(x*) + 4n l [ dx# - 2 f dx{x* 

Jo Jo 

x °(x*) = X* = x o(x*) ■ 


-x) 


d^ + ni 


(CIO) 


(Cll) 

(C12) 









We introduce the two transverse vectors e a = eg or e v which are orthogonal to the observed direction n Q , 


(C13) 

(C14) 

(C15) 


eg = (cos0cosy;,cos0sin(/?, —sin0), 
e v = (-siny;, cosy?, 0), 
n D = (sin 9 cos y>, sin 9 sin ip, cos 9). 

With this the deflection vector 5x l can be split into its radial part 5x r = n 0 i8x l and its transverse part 5x a = e°Sx'. 
From (1C11I) . we see that the transverse part contains one spatial gradients and consequently dominates over the radial 
part which has no gradient since the radial part of last term in (1C 1 111 is a total derivative. For the transverse part 
and the deflection angle we obtain to first order 

Sxi = -2 r dx(x*-x)e ai d^(x) 

Jo 

a? = -2 [ X ' d X ^^e ai dMx) • 

Jo X* 

3. Second order 

At second order the geodesic equation is 

d 2 x i _ _d^_k/_ _ r 
dx 2 dx (fc 0 ) 2 (fc 0 ) 2 

The first term on the right-hand side is negligible since it has less transverse derivatives than the second term. 
Indeed, at first order dk°/dx has no spatial derivative (see (1C9I) 1. and k l as only one transverse derivative (see (1C8[I ). 
In addition, when dk°/dx is taken at second order, it is multiplied by n l Q and therefore does not contribute to the 
transverse deviation Sx a . 

Neglecting the first term, we obtain 


(C16) 

(C17) 


(CIS) 


x l (x*) 


x b(x*) 



dx (fc °) 2 


x l o (x*) 



dx(x~X*) 


F 

(fc 0 ) 2 


(C19) 


where we neglect the perturbations at the observer. In (1C19I) . the integral is still performed on the perturbed geodesic. 
We rewrite this integral on the background geodesic by Taylor expanding the integrand around the background 
position. We get 


&4(x*) = 2 


d X ( X 

Jo 


X*) 


/r^r\ (2) 

V ( fc0 ) 2 ) 


6x( dj 


fv^k^\ (1) 
\ (fc °) 2 j 


Let us look at the first term 


/rpc^\ 

V (fc °) 2 ) 


( 2 ) 

— r* h v -u r i 

— 1 2nv^0 ^0 ' 1 1 jiv 



(C20) 


(C21) 


The second-order Christoffel symbols contain terms with one spatial derivative of the second-order gravitational 
potential. These terms are taken into account in standard lensing analyses, by using the halo-fit power spectrum in 
the first order expression for 5x a . In addition, the second-order Christoffel symbols contain coupling terms of the 
form i/jdiil}, with at most one spatial derivative. We can therefore neglect them. The second term in (IC21I) has at 
most two transverse derivatives, one in the first-order Christoffel symbols and one in the first order k\. It is therefore 
also negligible. The only relevant contribution comes therefore from the correction to the Born approximation, i.e. 
the second term in (IC20I) . We have 






















12 


Neglecting the second term we obtain for the deflection 

&4(x*) = - 4 [ dx(x* ^ x)dx{d 3 d l ^{x) • (C23) 

Jo 

The transverse part of Sx 2 (x*) 1S given by 

5x%(x*) =eiSx z 2 (x*) =8ef [ dx{x* - X)djd l ^/{x) [ dx'(x ~ xO&'&ix') ■ (C24) 

Jo Jo 

The derivative di can be expanded on the basis (n 0 , e#, e^) 

d J = n 0 o n k o d k + e J a e k d k . (C25) 

With this 

d : ^(x)dj^{x') = n k 0 d k ^{x)n l 0 dirf>(x!) + e k d k ^(x)e be de^(x) ■ (C26) 

Neglecting the first term, which has only radial derivatives, the transverse part of the second-order deviation vector, 
< 5 x 2 , and the second-order deflection angle, a 2 , become 

te%(x*) = 8 f dx f c^x'(x*-x)(x-x'K4'e ^ ’ J [9fca ^ 4'(x)]5 J N'(x , ) (C27) 

Jo Jo 

« 2 (x*) = 8 f^d X f X dx' ix * ~ x)i j ~ *' ) ete$e b i[d k d i *(x)]d j *(x') ■ (C28) 

Jo Jo X 

The radial part Sx^ = riddx^ has, by construction, less transverse derivatives and is therefore negligible with respect 
to the transverse part 6x?i ■ 


Appendix D: Details of the calculation of the lensed CMB Ci 
1 . a,urn in the full-sky 

We now use the deviation vector to calculate the full-sky di m at second order. The contribution from © reads 

/ /*x* rx 

dtt no Ye m (n 0 ) dx d X '{x*-x){x-x')e al eld i d j ^(n 0 ,x)e bk d k 'S(n 0 , X , )e^d m T{n 0 , x*) ■ (Dl) 
Jo Jo 


Using that 


ni mi bk ^ (A rn , mA\ti k , j k \ 

e e a e b e = -j;( e + e - +e+e_){e J + e_ + e J + e_), 


and the relations (1A17I) , (1A18I) and (1A19I) we obtain 

rX * f\ ,(x*-x)(x-xO 


/ rx* rx 

dfl no Y im (n 0 ) / dx dx' 

Jo Jo 


x 2 x'x* 


^ 2 4'(n Q , x) 0'&(n o , x') 0T{ n 0 , x*) 


(D2) 


(D3) 


+ 00^(n o ,x)0^(n o ,x')0T(n o ,x*) + $^(n 0 ,x)$^( n o,X / )$ T ( n o,X*)+ 0 ^( n o, x) #*(n 0 , x') 0T(n o , x*) . 

To calculate the angular power spectrum, we need to cross-correlate 8 2 a,i m with the unlensed ae> m ' 

ae'm'i n') = J dSl a > o Y e > m '(ri 0 )T(ri 0 , x*). (D4) 

Looking first at the first term in the square bracket of dD3l) we obtain 


(S 2 a£ m (n 0 )ae' m' «)} = 2 J dn n Y e m ) J dYl n iYf i rn i (n'jdx J dx 


x , ,(x*-x)(x-x') 


x 2 x'x* 


(D5) 


^ 2 ^(n 0 ,x) 0^{n o ,x')) ( 0T(n o ,x*) T ( n o,X*) 


X 
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where we have neglected the cross-correlation between the temperature and the gravitational potential coining from 
the integrated Sachs-Wolfe. The two-point function of the potential can be calculated using DU and the properties 
of the operators $ and $ on the Yi m ( n 0 ). We get 

^ 2 T(n 0 ,x)^(n 0 ,x')) = ~~ Y2 ^ ^ 2)l(f - l jl / ~f T *( k ’X)T*( k ,x')ji( k X)M k x') (D 6 ) 

x y^(-l) TO 2 ^m(n 0 )-iY^ m (n 0 ). 

m 

The sum over m can be simplified using 

^(- 1 ) ra 2 ^(n 0 )-ib, m (n„) = -^ 2 ^ m (n 0 ) 1 y/ m (n 0 ) = - 0) = 0 W. (D7) 

m m 

The same kind of calculation applies to the other terms in m leading to 

(5 2 a errL (n 0 )atm'(n' 0 )) = 0 ■ ( D8 ) 


2 . Ci in the flat-sky approximation 


Let us also briefly present the calculation of the lensed Ci in the flat-sky approximation, where the formalism is 
much simpler. As we will see, in the flat-sky approximation we find that the diagonal part of the two-point function 
vanishes in average. The off-diagonal part on the other hand vanishes under Limber approximation. But since in the 
flat-sky approximation the Limber approximation becomes exact, this shows that the flat-sky result is consistent with 
the full-sky result above. 

In the flat-sky approximation, the ai m are replaced by 2D multipoles 


m = 


d 2 yj c 

2-7T 


f(x 0 )e 


— if,-: 


(D9) 


where x G is the 2-dimensional observed position on the last scattering surface and t is the 2D variable of its Fourier 
transform. We want to calculate the contribution from the transverse part of S 2 T in (0 (remember that the radial 
part is negligible) 


S 2 T(£) = 


d 2 x c 

2 tt 


5x a d a T(x 0 )e 


-i£-z 


(DIO) 


Using (IC27I) for the second-order deflection vector and Fourier transforming (in 2D) the gravitational potential and 
the temperature, we obtain 

" d 2 x„ 


5 2 T{£) = 8 


° e -i£-x 0 


>0 


2 tt 

r+4i)-x, 


x 'd X 

Jo X*XX 


2n 


d 2 t 2 

27T 


dH 3 

27T 


(*1 -Wl •^(< 1 , X )^ 2 , X , )^ 3 , X .) 1 


(Dll) 


where we have used that in the flat-sky approximation the positions at which the potential is evaluated, x and x', are 
related to the observed position x D through 


X A I X 

x = —x„ and x = —x Q . 


X* 


The integral over x Q then gives 


d 2 x, 


2 t r 


oi(i 3-/+« i *+e 2 ^)-*o = 


X* 


2-K8[l z -l + h—+l 2 — 


X* 


X* 


(D12) 


(D13) 


With this (ID11I) becomes 


S 2 T(£) = 8 r d x [ 
Jo Jo 


dX 


,(x* - x)(x- x') f X* 


X*XX 


X 


2t r 


27T 


£-£s~£ 2 — I • t 

X* 


i-lo-i-A — 

X* 




^-£3-^)|,x)^, X')T(£ 3 , X *) • 


(D14) 
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We then cross-correlate S 2 T(£) with the unlensed multipole T(£') 

x * , f x , ,(x*-x)(x-x') f x*\ 3 [ d 2 £ 2 [ d 2 £ 3 


(6 2 T(£)T*(£')) = 8 f * d X f d X ‘ 
Jo Jo 


X*XX 


X 


2n / 2ir 


£-£ 3 -£ 2 — )■£ 

X* 


Using that 


<: T{£ 3 )T*{£')) = C T (£')6(£ 3 - (!) (! = \C \, 


and 


4- 


we obtain for the two-point correlation function 


(S 2 T(£)T*(£')) = 8 C T (C) [ * d X [ doc 

Jo Jo 


i\ [ x * ,, f x , ,(x*-x)(x-xO f x* 


x*xx 


X 


d 2 £ 2 

(2n)i 


X* 


£-£'- £> — \-£ 2 [£-£'- £ 2 — ) ■ (! 


(D15) 


(t-la- ^ 2 ^) • £3 (T(£s)T*(£')) (*(ji-£ 3 - £ 2 ^^-, x '\ *(£ 2 , x ')^ 


(D16) 


^£-£ 3 -£ 2 ^X±,X^) *(£ 2 , X ')j=C*(( 2 ,x,x')d(£ 2 +^£-£ 3 -£a^'j j) ^ = |€ 2 |,, (D17) 


X* 


c*(£ 2 ,x,x')6 l £ 2 *—*- + (£-S)— ) • 


. X* 


X 


(D18) 


The diagonal part of the two-point function vanishes. Indeed when £ = (!, the Dirac-delta function forces either 
X = x'i i R which case (ID 1811 vanishes due to the kernel x ~ x! or ^2 = 0, in which case (ID18I) also vanishes, due to the 

factor [£ — £'— £ 2 ■ £ 2 - 

The non-diagonal part £ ^ £' vanishes due to Limber approximation. The angular power spectrum of the gravita¬ 
tional potential is indeed given by 

/ dh 

—Pin{k)T^(k,x)T^(k,x')je 2 (kx)je 2 (kx') , (D19) 

where P ln is the initial power spectrum and Tj< is the transfer function. In Limber approximation the integral over k 
is simplified using 


J dkk 2 f{k)ji{kx)jt{kx) - ^2<5(x - x)f , (D20) 

and the non-diagonal contribution vanishes due to the factor x ~~ x! ■ The flat-sky result is therefore consistent with 
the full-sky result. 
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